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If} ' Abstract. Most calculations of quantum corrections in supersymmetric theories are made 

£N| , with the dimensional reduction, which is a modification of the dimensional regularization. 

However, it is well known that the dimensional reduction is not self-consistent. A consistent 
regularization, which does not break the supersymmetry, is the higher covariant derivative 
regularization. However, the integrals obtained with this regularization can not be usually 
calculated analytically. We discuss application of this regularization to the calculations in 
supersymmetric theories. In particular, it is demonstrated that integrals defining the /3-function 
(-h are possibly integrals of total derivatives. This feature allows to explain the origin of the 

exact NSVZ /3-function, relating the /3-function with the anomalous dimensions of the matter 
superfields. However, integrals for the anomalous dimension should be calculated numerically. 
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l/~j . 1. Introduction 

In order to deal with divergent expressions in the quantum field theory, it is necessary 
to regularize a theory. A proper choice of a regularization can simplify the calculations 
or reveal some features of quantum corrections. Most calculations in the quantum field 
theory were made with the dimensional regularization [1] in M£-scheme [2]. However, the 
S^ dimensional regularization breaks the supersymmetry and is not convenient for calculations in 

supersymmetric theories. That is why most calculations in supersymmetric theories were made 
with the dimensional reduction [3]. For example, the /3-function in supersymmetric theories was 
calculated up to the four- loop approximation [U El El El [B] . After a special redefinition of the 
coupling constant [21 [TO] the result coincides with the exact NSVZ /3-function, proposed in Refs. 

nn na aa H3Q . 

However, it is well known that the dimensional reduction is not self-consistent [TOJ. Removing 
the inconsistencies one breaks the explicit supersymmetry [16\ I17| . Then the supersymmetry 
can be broken by quantum corrections in higher loops [1S\ [TOl |2TJ] . In the N = 2 SYM theory 
this already occurs in the three-loop approximation [18[ [20] , while in the A = 4 SYM theory 
the supersymmetry is not broken even in the four-loop approximation [21] , Thus, a problem of 
regularization in supersymmetric theories is rather nontrivial [22J. 

For supersymmetric theories one can use the higher covariant derivative regularization, 
proposed by A.A.Slavnov [23[ 124]. Different versions of this regularization for supersymmetric 
theories, which do not break the supersymmetry, were proposed in |25[ 126] , Unlike the 
dimensional reduction, the higher covariant derivative regularization is consistent. However, it 



was not often applied to explicit calculations of quantum corrections, because the corresponding 
integrals have very complicated structure, and it is not easy to calculate them analytically, 
especially in higher loops. Moreover, some theoretical subtleties can raise nontrivial questions 
even in the simplest calculations [271 EH [29] ■ 

However, we argue that for supersymmetric theories this regularization has some very 
attractive features and can be used for the calculations. Namely, the integrals defining the 
/3- function are integrals of double total derivatives [301 EH [321 [33], and one of them can be 
calculated analytically. As a result, the /3-function is related with the anomalous dimension, 
producing the exact NSVZ /3-function without redefinition of the coupling constant. In this 
paper we demonstrate how this can proved in N = 1 SQED in all loops and for the general 
renormalizable N = 1 SYM in the two-loop approximation. 

2. Quantum corrections in N=l SQED, regularized by higher derivatives 

2.1. Higher derivative regularization 

The action of the massless N = 1 SQED in terms of superfields [341 [35] is written as 

S = -^Re f d 4 x d 2 9 W a C ab W b + \ [ d 4 x d 4 8 (> e 2V > + 0*e" 2y 0) . (1) 

The theory is regularized by adding the term with the higher derivatives 

5 A = ^Re J d 4 xd 2 9W a C ab (R(^) - l)\V b , (2) 

where -R(O) = 1 and R(oo) = oo. For example, it is possible to choose R = 1 + d 2n /A 2n . The 
gauge is fixed by adding 

% = - ^ / ** ^ {vrQ ^d 2 v + vrQd 2 d 2 v) . (3) 

Then the propagator will contain large degrees of the momentum in the denominator, and all 
loop diagrams beyond the one-loop approximation become convergent. The remaining one-loop 
diagrams are regularized by inserting the Pauli-Villars determinants [36J 

H( [ DtfDfae&y* 1 (4) 

into the generating functional, where 

Si = j I d A x d 4 9 (^e 2V >j + <^e~ 2y <£j) + (- f d 4 x d 4 6 M^^ + h.c.) (5) 

and Yl c i = 1) Yli c iM 2 = 0. It is important that the masses Mj are proportional to the 
parameter A. 

2.2. Three-loop f3 -function 

In order to find the /3-function we consider 

if = ~ J -^d 4 eV(-p)d 2 IL l/2 V(p)d- 1 (a,n/p) (6) 

and calculate 
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In the three- loop approximation the result can be written as (R k = R{k 2 / A 2 )) 
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(8) 



From this expression we see that the /3-function is defined by the integrals of (double) total 
derivatives [30|, I3T1 [32] . Such a structure allows to calculate one of the loop integrals analytically 
and, thus, to reduce a number of integrations. This can be done using the identity 
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where / is a nonsingular function which rapidly tends to at the infinity. It is equivalent to the 
identity 
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) = I ^(/M-/(o))=- T ^/(o). 
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(2vr) 4 q 2 dq 2 ' 

(This is a total derivative in the four-dimensional spherical coordinates.) Taking the integrals 
in Eq. ([5]) and comparing the result with the expression for the two-loop anomalous dimension 
of the matter superfield 
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in the considered approximation we obtain the exact NSVZ /3-function for the N = 1 SQED 
WT 



/3(ao) = ^(l-7(ao))+0(a§). (12) 



This equality of the well-defined integrals is obtained without any redefinitions of the coupling 
constant. Thus, with the higher derivative regularization the NSVZ scheme can be naturally 
defined. However, it is much more complicated problem to calculate analytically the integrals 
for the anomalous dimension. Possibly, in the lowest orders this can be done analytically, but it 
seems that in higher loops for this purpose one should use numerical methods. 

2.3. The exact result 

For the N = 1 SQED it is possible to demonstrate that the features discussed in the previous 
section take place in all loops. In particular, the integrals defining the /3- function are integrals 
of double total derivatives, and the /3-function coincides with the NSVZ expression without a 
redefinition of the coupling constant [37] . 

For this purpose, first, we make the substitution 

v -»• e a e a 9% = e A , (13) 

which allows to extract the function d~ l in Eq. (jSJ). Then we will try to present the sum of 
Feynman diagrams as integrals of total derivatives, which in the coordinate representation are 
given by 

Tr ( [x^ , Something] ) = 0. (14) 

We start with the expression for the part of the effective action corresponding to the two-point 
Green function of the gauge superfield |37] 

Ar^ } = ^-2i(Tr(VJo*)) 2 -2iTr(VJ *VJ *)-2iTr(V 2 Jo*))+terms with * + (PV), (15) 

where 

1 ~- 1 

* = l-(e 2V -l)D 2 D 2 /Wd 2 ' * ~ 1 - (e~ 2V - l)D 2 D 2 /16d 2 ^ 

encode sequences of vertexes and propagators on the matter line, and Jo = e 2V D 2 D 2 /16d 2 is 
the effective vertex. (PV) denotes contributions of the Pauli-Villars fields. The first term in Eq. 
(|15p is a sum of diagrams in which external lines are attached to different loops of the matter 
superfields. The second term is a sum of diagrams in which external lines are attached to a 
single line of the matter superfields. The last term is not transversal. The sum of such terms 
vanishes due to the Ward identities. 

After substitution (|13|) and some algebraic transformations [37] the first term in Eq. (|15p 
gives the contribution 

- 2i ^x(( Tr ( - 2W^,in(*) - in(*)] + fe c {Y)c d e d [yl, H*) - H*)\) + (pv)) 2 ). (17) 

Similarly, the second term in Eq. (|15p gives 
. d 



rf ,„A W 



Vu> 



(y**)*, ln(*) + ln(*) ) + (PV) - terms with a (5-function. (18) 



The third term in Eq. f)15f) vanishes after substitution (|13p . From expressions (|17p and f)18|) we 
see that in all orders the /3-function is given by integrals of double total derivatives. A different 
method to see this [31] is based on the covariant Feynman rules in the background field method 



Expressions (|17p and (|18p can be calculated explicitly in the three-loop approximation. The 
result coincides with Eq. (JSj) - 

Terms with the 5-function in Eq. (|18|) appear due to the identity 

W, %\ = H^, -|d = -2^ 2 ^ 4 fe) = -2vr 2 ^(p). (19) 

Due to this <5-function one of loop integrals can be calculated and a number of integrations is 
reduced. Qualitatively, we consider all diagrams in which two external gauge lines are attached to 
the same graph. Integration of the <5-function corresponds to cutting a matter line in this graph 
|31| . This gives diagrams with two external matter lines, defining the anomalous dimension. 
Thus, the /3-function in a certain loop is reduced to the anomalous dimension in the previous 
loop. These arguments can formulated rigorously [37] in all orders and allow to obtain the exact 
NSVZ /3-function 

/3(a) = ^(l - 7 (a)). (20) 

Note that deriving Eq. (|20p one does not redefine the coupling constant, as in the case of 
the dimensional reduction [9j[T0]. Therefore, with the higher derivative regularization we can 
naturally define the NSVZ scheme for the N = 1 SQED. 

3. Two-loop /3-function with the higher covariant derivative regularization in the 
non-Abelian case 

Let us consider a general renormalizable iV = 1 super symmetric Yang-Mills theory with a 
gauge group G and matter superfields <j>i in a representation R, in the massless limit: 

S= -^Retr J d 4 xd 2 9W a C ab W b + ^ f ' d 4 xd A 9 (0*)*(e 2V V >i + 

+ (^ J d 4 xd 2 9\ ijk (t>i(f> j (j> k +h.cA (21) 

This theory is invariant under the gauge transformation if 

{T A ) m i X mjk + (T A )J\ imk + (T A ) m k X ijm = 0. (22) 

Below we assume that this condition is satisfied. In order to introduce the higher covariant 
derivative regularization and calculate the /3-function in this case we will use the background 
field method. In the supersymmetric case [34|, 135] we split the gauge superfield (which is 
below denoted by V) into the quantum part V and the background field ft according to the 
prescription e = e e e . Then the gauge can be fixed without breaking the background 
gauge invar iance: 

5 Ef = 3_tr fd A xd 4 e(VD 2 D 2 V + VD 2 D 2 v). (23) 

32e^ J V / 

(In our notation D, D, and D a are the background covariant derivatives.) Certainly, the 
gauge fixing procedure also requires introducing the Faddeev-Popov and Nielsen-Kallosh ghosts. 



The higher covariant derivative regularization can be also introduced without breaking the 
background gauge invariance. This can be done by different ways. For example, it is possible to 
add 



S A = — tr Re / d A x d A 6 V^-%1 V + - [ d*x d 
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to the action, assuming that n > m. (It is important that for a theory with a nontrivial cubic 
superpotential a term with the higher covariant derivatives should be also introduced for the 
matter superfields.) 

As in the case of N = 1 SQED, the higher covariant derivative term does not remove 
divergences in the one-loop approximation, and in order to regularize them the Pauli-Villars 
determinants should be inserted into the generating functional. The Pauli-Villars fields should 
be introduced for the matter superfields and all ghosts. (A contribution of the gauge superfields 
in the one- loop approximation vanishes.) Masses of the Pauli-Villars superfields (f>i are 
proportional to the parameter A: Mj J = af A and satisfy the relation M l j (Mj)jk = Mj5 k . 

Using the background gauge invariance it is possible to choose fi = fi + = V. Then the 
function d" 1 is defined by 

r S/ 2) =-^tv J 0^d^V(-p)d 2 U 1/2 V(p)d- 1 (a,\^/p). (25) 

In order to find the f3- function we again use prescription ([7D and find the derivative of d _1 with 
respect to In A in the limit of the vanishing external momentum. In the non-Abelian case the 
result is given by 
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f ACD f BCD^ C2d AB. r ^ 5 AA . ' (27) 

The integrals defining the /3-function are given by 

li = I 4 (0) - 5>jii(Mi), i = 0,2,3 ? (28) 

i 

where for simplicity we do not write the ghost constributions (they are also given by integrals 
of double total derivatives) and 
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Thus, we see that the /3-function is given by integrals of double total derivatives. In the 
considered (two- loop) approximation they can be easily calculated analytically using Eq. ([9]): 



- , „ , a 3 / _ 2 , m ,„^ , 2„, n , <„, n , A a 2 C(R)^X* kl X ikl 



/3 = - - (3C 2 - T(RJ) + ^_ ( - 3Cf + T(i?)C 2 + -C(UVC(«) j , 8vr 3 r 

(34) 
This expression should be compared with the one-loop anomalous dimension 

y (o) __2^ + ^! + ... (35 , 

Then we see that in the considered approximation the /3-function agrees with the exact NSVZ 
/3-function 

a 2 [3C 2 - T(R) + C(2?y 7 /(«)A 

«°> = mi- W "• (36) 

4. Conclusion 

Although it is generally believed that the integrals appearing with the higher covariant 
derivative regularization have too complicated structure, we see that in the super symmetric 
case some of them can be calculated analytically. This makes possible analytical multiloop 
calculations in supersymmetric theories with this regularization. In principle, in the lowest 
loops it is not very difficult to construct integrals corresponding to various Green functions. 



A very attractive feature of the higher covariant derivative regularization is that all integrals 
defining the /3-function in the super symmetric case seem to be integrals of double total 
derivatives. As a consequnce, one of them can be calculated analytically. (For N = 1 SQED 
this was proved exactly in all loops. For the general renormalizable N = 1 SYM theory this was 
so far verified only in the two- loop approximation.) In both considered cases the factorization 
of integrands into total derivatives allows to obtain the exact NSVZ /3-function, and for N = 1 
SQED for this purpose it is not necessary to redefine the coupling constant. (In the non-Abelian 
case we cannot so far make this conclusion, because the calculation was made only in the two- 
loop approximation, where the /3-function is scheme- independent.) 

Acknowledgments 

This work was supported by Russian Foundation for Basic Research grants No 11-01-00296-a 
and 11-02-08451-z. I would like to thank the organizers of the ACAT conference for supporting 
my participation. I am also very grateful to prof. A.L.Kataev for valuable discussions. 

References 



[l] 
[2] 
[3] 
[4] 

[5] 
[0] 
[7 

[s; 

[io; 

[ii 

[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20' 
[21 
[22 
[23" 
[24 
[25" 
[26 
[27 
[28 
[29 
[30" 
[31 
[32 
[33 

[34; 

[35 

[36" 
[37 
[38 
[39" 



t'Hooft G and Veltman M 1972 Nucl.Phys. B 44 189. 

Bardeen W A, Buras A J, Duke D W and Muta T 1978, Phys.Rev. D18 3998. 

Siegel W 1979 Phys.Lett. B 84 193. 

Ferrara S and Zumino B 1974 Nucl.Phys. B 79 413. 

Jones D R T 1975 Nucl.Phys. B 87 127. 

Avdeev L V and Tarasov O V 1982 Phys.Lett. B 112 356. 

Jack I, Jones D R T and North C G 1996 Phys.Lett B 386 138. 

Harlander R V, Jones D R T, Kant P, Mihaila L and Steinhauser M 2006 JHEP 0612 024. 

Jack I, Jones D R T and North C G 1997 Nucl.Phys. B 486 479. 

Jack I, Jones D R T and Pickering A 1998 Phys.Lett. B 435 61. 

Novikov V A, Shifman M A, Vainshtein A I and Zakharov V I 1983 Nucl.Phys. B 229 381. 

Novikov V A, Shifman M A, Vainshtein A I and Zakharov V I 1985 Phys.Lett.^, 166 329. 

Shifman M A and Vainshtein A I 1986 Nucl.Phys. B 277 456. 

Vainshtein A I, Zakharov V I and Shifman M A 1985 JETP Lett. 42 224. 

Siegel W 1980 Phys.Lett. B 94 37. 

Avdeev L V, Chochia G A and Vladimirov A A 1981 Phys.Lett. B 105 272. 

Stockinger D 2005 JHEP 0503 076. 

Avdeev L V 1982 Phys.Lett. B 117 317. 

Avdeev L V and Vladimirov A A 1983 Nucl.Phys. B 219 262. 

Velizhanin V N 2009 Nucl.Phys. B 818 95. 

Velizhanin V N 2011 Phys.Lett. B 696 560. 



Jack I and Jones D R T 1997 Regularisation of supersymmetric theories Preprint hep-ph/9707278 

Slavnov A A 1971 Nucl.Phys. B 31 301. 

Slavnov A A 1972 Theor.Math.Phys. 13 1064. 

Krivoshchekov V K 1978 Theor.Math.Phys. 36 745. 

West P 1986 Nucl.Phys. B 268 113. 

Martin C and Ruiz Ruiz F 1995 Nucl.Phys. B 436 645. 

Asorey M and Falceto F 1996 Phys.Rev. D 54 5290. 

Bakeyev T and Slavnov A A 1996 Mod.Phys.Lett. A 11 1539. 

Soloshenko A A and Stepanyantz K V 2004 Theor.Math.Phys. 140 1264. 

Smilga A and Vainshtein A 2005 Nucl.Phys. B 704 445. 

Stepanyantz K V 2005 Theor.Math.Phys. 142 29. 

Pimenov A B, Shevtsova E S and Stepanyantz K V 2010 Phys.Lett. B 686 293. 

West P 1986 Introduction to super-symmetry and supergravity (World Scientific). 

Buchbinder I L and Kuzenko S M 1998 Ideas and methods of supersymmetry and supergravity (Bristol and 

Philadelphia, Institute of Physics Publishing). 
Faddeev L D and Slavnov A A 1990 Gauge fields, introduction to quantum theory (Benjamin Reading). 
Stepanyantz K V 2011 Nucl.Phys. B 852 71. 
Grisaru M T and Zanon D 1985 Nucl.Phys. B 252 578. 
Grisaru M T, Milewski B and Zanon D 1986 Nucl.Phys. B 266 589. 



